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ABSTRACT 

Using current models of spectrophotometric properties of single age, single metallicity 
stellar populations I have computed the Mg2 , H/3 , Fe52 and Fe53 line strengths for 
stellar populations with a metallicity spread. The comparison of these models with 
the nuclear indices of early type galaxies yield the following major conclusions. The 
metallicity distribution of the closed box, simple model for the chemical evolution of 
galaxies is not able to account for Mg2 and Fe52 , Fe53 values in excess of ^ 0.27,3 and 
2.7, respectively, which are observed in the nuclei of a large fraction of Ellipticals. To 
reproduce the line strengths in these galaxies high average metallicities, small metallicity 
dispersion and old ages are required. In particular, Mg2 values of ~0.3 are reproduced 
only with a metallicity distribution ranging from O.5Z0to ~ 3Zq, and ~ 15 Gyr old 
stellar populations. I interpret the data as indicating that the gas out of which the 
nuclei of ellipticals formed was pre-enriched, to larger metallicities for increasing Mg2 . 
The presence of a metallicity dispersion does not alter the relation between Mg2 and 
Iron indices with respect to the SSP models. Thus, the need for a Mg/Fe overabundance 
in the strongest lined galaxies is confirmed, and I present a simple way to estimate the 
[Mg/Fe] ratio on the basis of existing models with solar abundance ratios. 

Key v^rords: Galaxies: elliptical & lenticular, stellar content, metallicity, formation 



1 INTRODUCTION 

The question of the age of the bulk of the stars in ellipti- 
cal galaxies is still subject to a large debate. Opposed to 
the classical picture in which Ellipticals are basically inhab- 
ited by ~ 15 Gyr old stellar populations (see e. g. Renzini 
1986), O'Connell (1986), amongst others, proposed that a 
substantial component of stars as young as 5 Gyr has to be 
present to account for the observed spectral energy distribu- 
tion of Ellipticals. In more recent years many observational 
evidences have been found to support the notion of elliptical 
galaxies formation at high redshift, including the tightness 
of the colour— central velocity dispersion (a) relation found 
for Ellipticals in Virgo and Coma by Bower, Lucey & Ellis 
(1992) ; the thinness of the fundamental plane (Renzini & 
Ciotti 1993) for the Ellipticals in the same two clusters; the 
modest passive evolution measured for cluster Ellipticals at 
intermediate redshifts (Franx & van Dokkum 1996, Bender, 
Ziegler & Bruzual 1996); the negligible luminosity evolution 
observed for the red galaxies (Lilly et al. 1995) and for early 
type galaxies in clusters in the redshift range z < 1 (Dickin- 
son 1996); the detection of bright red galaxies at redshifts as 
large as 1.2 (Dickinson 1996). On the other hand, hints for a 
continuous formation of Ellipticals in a wide redshift range 



have also been found: the relatively large H/3 values mea- 
sured in a sample of nearby Ellipticals, which could indicate 
a prolongued star formation activity in these galaxies, up to 
~ 2 Gyr ago (Gonzalez 1993, Faber et al. 1995); the appar- 
ent paucity of high luminosity Ellipticals at z~ 1 compared 
to now (Kauffmann, Chariot & White 1996). 

Two competing scenarios have been proposed also for 
the process leading to the formation of the bulk of the stars 
in Ellipticals: early merging of lumps containing gas and 
stars (e.g. Bender, Burstein & Faber 1993), in which some 
dissipation plays a role in establishing the chemical structure 
of the outcoming galaxy; and the merging of early formed 
stellar systems, occurring in a wide redshift range, and pref- 
erentially at late epochs, following the hierarchical formation 
of structures (Kauffmann, White & Guiderdoni 1993). 

In order to help understanding when and how ellipti- 
cal galaxies formed I have computed synthetic spectral in- 
dices for stellar populations with a metallicity spread, and 
compared them to the corresponding observations in the nu- 
clei of Ellipticals. The study of line strengths in the spec- 
tra of early type galaxies has shown to be a powerful tool 
for investigating the age and the metallicity of these sys- 
tems (Faber et al. 1995; Fisher, Franx and Illingworth 1995; 
Buzzoni 1995b and references therein) . With few exceptions 
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(Vazdekis et al. 1996, Bressan, Chiosi & Tantalo 1996), most 
of the authors have interpreted the observed Hne strengths 
through comparisons with theoretical models constructed 
for single age, single metallicity stellar populations (SSPs). 
The major results of these studies can be sumarized as fol- 
lows: 

i) the Mg2 indices measured in elliptical galaxies are con- 
sitent with the notion that these systems are inhabited by 
old stellar populations, the differences in Mg2 tracing dif- 
ferences in average metallicity (Buzzoni, Gariboldi & Man- 
tegazza 1992). However, the difficulty in determining sepa- 
rately age and metallicity (Renzini 1986) weakens consider- 
ably this simple picture (Worthey 1994); 

ii) the H/3 line strength offers an opportunity to break the 
age— metallicity degeneracy, if this index is measuring the 
temperature of turn-off stars (see Faber et al. 1995). In 
this view, the data derived for a sample of nearby Ellipti- 
cal galaxies indicate that the ages of the stellar populations 
in their nuclei span a wide range, the weakest Mg2 galaxies 
being the youngest (Gonzalez 1993, Faber et al. 1995) ; 

iii) the Magnesium to Iron abundance ratio in the nuclei of 
the highest Mg2 Ellipticals is likely to be larger than solar 
(Gorgas, Efstathiou & Aragon Salamanca 1990; Worthey, 
Faber & Gonzalez (1992), hereinafter WFG; Davies, Sadler 
& Peletier 1993). Weiss, Peletier and Matteucci (1995) esti- 
mate [Mg/Fe] ranging from 0.3 to 0.7 dex within the bright- 
est ellipticals. 

Real galaxies host composite stellar populations, with 
a spread in the major parameters like age and metallicity. 
This may apply also to the nuclei of galaxies, where typically 
~10^ LQare sampled, corresponding to ~100 bright globular 
clusters. The existence of a substantial metallicity spread in 
elliptical galaxies and bulges is supported by direct observa- 
tions. For example, the Colour-Magnitude diagram of a field 
in M32 (Freedman 1989) shows a wide red giant branch, 
corresponding to stars spanning a metallicity range of 0.6 
dex approximately. The K-giants in the galactic bulge have 
metallicities ranging from ^ O.IZqIo ~ 5Z0(Rich 1988). Fi- 
nally, the mere evidence for abundance gradients in elliptical 
galaxies, as inferred from line strengths gradients, indicates 
that a metallicity spread is present in these systems. This 
last argument applies to galaxies as a whole: whether or not 
their nuclei host stellar populations with a spread in metal 
content depends on the modalities of the galaxy formation. 
However, due to projection effects, a substantial fraction of 
the light measured in the nuclei of Ellipticals comes from re- 
gions located outside the three dimensional core. This frac- 
tion can e.g. amount to ~ 50 % for King models (Binney & 
Tremaine 1987) . Therefore a pure radial metallicity gradient 
translates into a metallicity spread in the stellar population 
contributing to the light measured in the galactic nuclei. 

In this paper, the effect of a metallicity spread on the 
integrated indices is investigated, viewing a given galaxy (or 
a portion of it) as the sum of SSPs. These models are com- 
pared to the relevant observations to derive conclusions on 
the stellar content of the nuclear regions of early type galax- 
ies and inferences on their formation process. Section 2 de- 
scribes how the models are computed, and the results are 
presented in Section 3. In Section 4 the models are com- 
pared to the observational data, and in Section 5 the impli- 
cations for the formation of elliptical galaxies are discussed. 
The main conclusions are sumarized in Section 6. 



2 COMPUTATIONAL PROCEDURE 

The spectral indices considered in the present work are de- 
fined as follows (Burstein et al. 1984): 
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where the various Fi denote the fluxes measured within 
the spectral windows of the different lines, centered at 
A ~ 5175, 4863, 5267 and 5334 A for the Mg2 , H/3 , 
Fe52 and Fe53 features, respectively. Fc are the pseudocon- 
tinuum fluxes measured at the line location, as interpolated 
from the fluxes measured in adjacent windows, and Ap, 
Apo are the wavelength widths of the windows in which the 
H/3 ad Iron indices are measured. 

The above definitions apply to the spectra of single 
stars, of SSPs and of collections of SSPs, inserting the ap- 
propriate values for the fluxes. Therefore, for a collection 
of N simple stellar populations, each contributing a fraction 
<I>s to the total bolometric flux Fboi of the composite stellar 
population, the integrated indices are given by equations (|l|) 
to (0) with 
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where the subscript S refers to the single SSP. The spectral 
energy distribution of an SSP, and particularly the various 
flux ratios, are controlled by a number of parameters, in- 
cluding the metallicity (Z), age (t), helium content {Y) and 
the elemental abundances of the population. Thus, the in- 
tegrated indices of a collection of SSPs depend on how the 
fractionary bolometric flux $s is distributed over the range 
covered by all the relevant parameters. The problem of de- 
riving these parameters from the observed line strengths can 
be symplified considering various suitable indices, each con- 
trolled by different parametes. For example, Gonzalez (1993) 
uses a combination of Magnesium and Iron line strengths 
(mostly sensitive to Z) and the H/3 index (mostly sensitive 
to t) to determine age and metallicity of a sample of early 
type galaxies. Still, in order to map the integrated indices 
of collections of SSPs into their fundamental properties one 
needs to account for the presence of a possible spread in the 
parameters which control the various line strengths. 

One possible approach to this problem consists in com- 
puting models for the chemical evolution of galaxies, which 
automatically yield the distribution of the SSPs over the 
fundamental parameters (e.g. Vazdekis et al. 1996, Tantalo 
et al. 1996). The output of these models, though, depends 
on the specific ingredients used, like the adopted star for- 
mation rate, initial mass function, ratio of dark to luminous 
matter, nucleosynthesis, criteria for the establishment of a 
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galactic wind, etc. A different approach consists in explor- 
ing the dependence of the various indices on the presence 
of a spread in the populations parameters adopting a phys- 
ically motivated function <E>s, and using relations (hi) to M). 
This approach, which has the advantage of allowing an easy 
exploration of the parameter space by simply changing the 
$s functions, will be adopted here. Also, I will restrict to 
considering collections of SSPs all with the same age, but a 
substantial spread in metallicity. The results are then meant 
to describe the effects of the presence of a metallicity spread 
in a stellar population formed within a short time scale, so 
that the integrated indices are not appreciably inffuenced by 
age differences in the individual components. 

2.1 Line strengths for simple stellar 
populations 

In order to compute the integrated indices of composite stel- 
lar populations one has to know the Fi/Fhoi and Fc/Fhoi ra- 
tios of the SSPs as functions of Z and t. Available SSP mod- 
els in the literature tabulate bolometric corrections, colours 
and spectral indices for different ages and metallicities. I 
then write: 



^ Ls ^ _L Ms Mr 
^ ^ Lt ^ M Mt Lt' 



(8) 



-Fbol ^ Fc^ ^ Fbol ^ 



(7) 



for each SSP, and approximate Fc with the flux in the V 
band, for Mg2 , Fe52 and Fe53 , and with the B band flux 
for H/3 . Although the pseudocontinuum fluxes do not corre- 
spond precisely to the flux in the V or the B band, the main 
conclusions of this paper are not affected by this approxi- 
mation. For example, the integrated Mg2 indices computed 
with Fc = Fb differ from those computed with Fc — Fv by 
less than 2 percent. 

Two sets of SSP models are used here: Buzzoni's mod- 
els (Buzzoni et al. 1992; Buzzoni, Mantegazza & Gariboldi, 
1994; Buzzoni 1995a), and Worthey's models (Worthey 
1994), hereinafter referred to as B and W respectively. The 
metallicity range covered by W models goes from O.OlZQto 
~3Z0 (Z0~ 0.017), while B models encompass a larger Z 
range, from ~ O.OOGZqIo ~ 6Zq. Each metallicity is char- 
acterized by one value for the helium abundance: Worthey 
(1994) assumes Y = 0.228 -I- 2.7Z, while in B models Y in- 
creases less steeply with Z, according to K ~ 0.23 -I- Z. Be- 
sides, the isochrones used to construct the SSPs have solar 
abundance ratios. This corresponds to specifying the chem- 
ical trajectory followed in the evolution of the composite 
stellar population. 

The two sets of models present systematic differences in 
the bolometric output, broad band colours and line strengths 
(Worthey 1994, Buzzoni 1995b). At least part of these differ- 
ences can be ascribed to the different choices of the AY / AZ 
parameter (Renzini 1995) and to the different fltting func- 
tions (i.e. the dependence of individual stellar indices on 
effective temperature, gravity and metallicity) adopted in 
the computations. 

2.2 The metallicity distribution functions 

The contribution of a given SSP to the bolometric light of 
a composite stellar population of total mass Mt and total 
luminosity Lt can be written as 



The distribution function $s is then proportional to the 
mass distribution over the metallicity range, through the in- 
verse of the {M/L)s ratio, which depends on the metallicity 
and helium content (Renzini 1995). 

The closed box, simple model for the chemical evolution 
of galaxies predicts that the distribution of the stellar mass 
over the total metallicity follows the relation / [Z) oc e~^^^ 
(Tinsley 1980), where y is the stellar yield. According to 
this model, most stars are formed at the lowest metallici- 
ties. More sophisticated models for the chemical evolution 
of elliptical galaxies, which take into account the occurrence 
of galactic winds (Arimoto & Yoshii, 1987; Matteucci & Tor- 
nambe, 1987), also predict the existence of a substantial 
metallicity spread in the stellar content of elliptical galax- 
ies. In these models, the more massive the galaxy, the later 
the galactic wind sets in, and further chemical enrichment 
is achieved. Correspondingly, as the galactic mass increases, 
the metallicity distributions get skewed towards higher Z 
values. Nevertheless, a substantial fraction of low Z stars is 
always present, and indeed Arimoto & Yoshii's metallicity 
distributions (by number) for model Ellipticals with mass 
ranging from 4 x 10^ to 10^^ Mq are well described by a 
closed box model relation : 



f{Z) 



exp(-Z/y) 
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with a minimum metallicity Zm — O.OIZq, a maximum 
metallicity Zm increasing from ~ 2Z0to ^ OZ©, and y vary- 
ing from 2ZQto SZq. So, the metallicity range spanned by 
the stars in these model ellipticals goes from ~ 2.3 to ~ 2.8 
dex. 

Relation (^ flnds also observational support from the 
direct determination of the metallicity distribution of K- 
giants in the bulge of our galaxy, which appear to fol- 
low a close box model relation with Zm ~ 5Z0and y ~ 
2ZQ(Rich 1988). Mc William and Rich (1994) revised Rich 
(1988) metallicity scale for [Fe/H] towards values lower by 
~ 0.3 dex, but flnd a [Mg/Fe] overabundance of the same 
amount, and state that the distribution of [(Fe-|-Mg)/II] in 
the bulge stars may agree with Rich (1988) [Fe/H] distribu- 
tion. 

I then adopt eq. to describing the Ms/Mt distri- 
bution over the metallicity for the composite stellar popula- 
tions, and explore the effect of different values for the three 
parameters Zm, and y. For a fixed (low) Z^, increas- 
ing values of Zm (and yield) are meant to describe stellar 
populations produced in an environment in which the chem- 
ical processing is terminated at progressively higher levels of 
completion. These sequences of models, then, conform to the 
predictions of chemical evolution models with galactic winds 
for increasing galactic mass. Different values of Zm, instead, 
characterize different degrees of pre-enrichment of the gas. 

To derive the $s distribution the behaviour of the 
(M/L) ratio for the SSPs with increasing metallicity and 
helium content needs to be specified. In B models, which are 
characterized by a low AY / AZ parameter, M/L increases 
with Z, going from 2.2 to 4.5 for Z ranging from O.OlZ0to 
3Zq, for the 15 Gyr old SSPs. In W 17 Gyr old models, 
the mass to (bolometric) luminosity ratio increases mildly 
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with metallicity up to a maximum value of 4.2 reached 
at O.5Z0,and it decreases afterwards, down to 3.3 at 3Zq. 
While the different behaviour reflects the different AY/AZ 
(Renzini 1995), these M/L values are not directly compa- 
rable, the total mass M being computed with different pre- 
scriptions (Worthey 1994). Besides this, neither of the two 
M/L correspond to what should be inserted in eq. (^) : Buz- 
zoni's M values do not take into account the mass locked 
into stellar remnants; Worthey's M values do not include 
the remnant masses for stars born with M >2 Mq, but do 
include both the remnant and the returned mass for stars 
born with mass between the turn-off mass and 2M0 . I then 
chose to neglect the dependence of (M/L)s on the metal- 
licity of the SSPs, and adopt "^s oc f{Z). If AI/L increases 
with metallicity, as in B models, this approximation leads to 
an overestimate of the contribution of the high metallicity 
SSPs in the integrated indices for the composite stellar pop- 
ulations. If, on the contrary, the M/L given in W models 
are more appropriate, the contribution of the high Z popu- 
lations will be underestimated. Notice, however, that in W 
models AI/L varies by only a factor of 1.2 for Z varying 
from O.OlZ0to 3Z0. In the following, composite stellar pop- 
ulations with this $3 distribution function will be briefly 
referred to as CSPs. 



3 RESULTS OF THE COMPUTATIONS 

Following the prescriptions described in Sect. 2, the inte- 
grated spectral indices are given by the following relations: 



Me 



H/3 = 



-2.5 Log 



-z/v 



J^^' dZ 10»'* ^C'- 



/I" dZ H/3 s 10° 



dZ 10° * e-z/y 



Fe52 = 



Fe53 = 



/I" dZ Fe52 s 10° 



/^^^' dZ 10°-* ^^v e-z/y 



/I" dZ Fe53 s 10° 



/I" dZ 10°-* ^^v e-z/y 



(10) 



(11) 



(12) 



(13) 



where BCb and BCv denote the bolometric corrections to 
the B and the V band magnitudes, respectively, and the S 
index denotes the SSP quantities, which depend on metal- 
licity and age. One can immediately notice that, since high 
Z populations yield a relatively low contribution in the B 
and V bands, the low metallicity component of a compos- 
ite stellar population will tend to dominate the integrated 
indices. 

The behaviour of these indices as functions of the av- 
erage metallicity of a composite stellar population has been 
explored by considering different values for Zm and Zm- In 
order to characterize each distribution in terms of one pa- 
rameter, an average metallicity has been computed, defined 
by the following relation: 



[< Fe/H >] = Log 



(Z/^)0 



-ziv 



J^^dZ e-^iy 
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Other definitions of average metallicity can be found in the 
literature (see e.g. Arimoto & Yoshii 1987). In this respect, 
it should be noticed that the quantity < [Fe/H] > differs 
sistematically from [<Fe/H>] given by eq. ( |l4| ) because the 
former assigns more weight to the low metallicity tail of 
the distribution. The difference between the two quantities 
amounts to ~ 0.2 dex for the widest Z distribution consid- 
ered here (O.OIZ0 <Z < 6Z0). The choice of using eq. (|l4|) 
is motivated by the fact that this parameter better describes 
the mass fraction of metals in the CSP, for the given /(^) 
distribution function. 



3.1 Indices versus metallicity using B models 

I will now describe the results of the integration of equations 
( p^ to (^3|) for the different Z distributions, using B and W 
sets of SSP models. Among the various models by Buzzoni 
I have considered those computed with Salpeter IMF and 
red horizontal branches. Analogous assumptions character- 
ize Worthey's models. 

Figure 1 displays the integrated indices obtained using 
B SSPs at 15 Gyr, which are shown as thick lines in the four 
panels. The thin lines illustrate the effect of the presence of a 
metallicity spread (see captions), the lowest of which corre- 
sponds to Zm = O.OIZ0, and shows the expected behaviour 
of the indices for galaxies of increasing mass, in the frame 
of the wind models for the chemical evolution of galaxies. 
The different lines shows the results obtained with differ- 
ent Zm, describing the effect of assuming different degrees 
of preenrichment of the gas. 

The loci described by the CSPs for increasing average 
metallicity are shallower than the pure SSP relations: this 
is due to the fact that the higher Z populations contribute 
less than the lower metallicity ones in the optical bands, 
where the considered spectral indices are measured. As a 
consequence, the Mg2 and Iron indices of these CSPs never 
reach values as high as those characteristic of the highest Z 
SSPs, unless the metallicity spread is extremely small. Had I 
used the B-band fiux as a measure of the pseudocontinuum 
for calculating the integrated Mg2 index the effect would 
be stronger, since the high metallicity populations would re- 
ceive even less weight. Taking into account the [M / L)s ratio 
dependence on Z as given in B models, when computing $s, 
would also strengthen the difference between SSP and CSP 
models. 

At any given average metallicity, the Mg2 and Iron in- 
dices for CSPs are weaker (and H/3 is stronger) than the 
corresponding values for SSPs. For example, a value of 
[<Fe/H>] = 0.22 can be obtained with a composite stellar 
populations with parameters (Zm,ZM,J/)= (O.Ol,5,2)xZ0. 
The differences between the line strenghts of such a com- 
posite population and those of the SSP with the same 
[Fe/H] amount to A Mg2 ~ -0.05, A H/3 ~ 0.18, A Fe52 ~ 
—0.5 and A Fe53 ~ -0.4. These results refer to the y = 
2xZ0case. Adopting y — 4 XZ0 (dotted lines), so as to en- 
hance the fraction of the high Z component, the differences 
are only slightly smaller. This effect is particularly impor- 
tant when dealing with the highest metallicity galaxies. 

In much the same way, at any given value of the 
considered spectral index, the average metallicity of 
a composite stellar population is higher than the 
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Figure 1. Model line strenghts using Buzzoni's 15 Gyr old SSPs, as functions of metallicity. The thick lines are the loci of pure SSP 
models. The thin lines display models for CSPs for various values of the Zm parameter. The following cases are shown: Zni/ZQ = 
(0.01,0.1,0.3,0.5,1,1.5,2), and the filled squares mark the values for SSPs with Z = Z^- Along each line, Zm varies from up to 0.1. 
The filled circles mark CSP models with Zm = O.IZq, and Zm = (O.5,l,1.5,2,4,6)xZ0. The effect of different values of the parameter y 
is also shown. Finally, the dot-dashed lines show the typical ranges spanned by observational data for the nuclei of elliptical galaxies. 



metallicity of the SSP which has the same index. This 
is illustrated in Figure 2, where I plot, as a function of Mg2 , 
the difference (A[Fe/H]) between the average metallicity of 
the CSPs and that of SSP models with the same value of 
the Mg2 index. Along any line, the metallicity distributions 
have the same Zm and increasing Zm, as in Figure 1. Figure 
2 shows that the metallicity inferred from a given Mg2 value 
using SSP models is lower than what would be derived using 
CSP models, the difference being larger the wider the metal- 
licity distribution. The lower weight received by the high Z 
populations in eq. (|l^) causes the rapid growth of A [Fo/H] 



as the high Mg2 ends of the curves are approached: high 
Mg2 hne strengths are obtained only with very large Zm 
values. This is not a small effect: at Mg2 = 0.26 the differ- 
ence in the metallicities amounts to 0.3 dex, for a metallicity 
distribution extending down to Zm ~ O.OIZ©. 

Adopting a large y the effect is still important, while it 
gets small only reducing the widths of the metallicity distri- 
butions, (i.e. increasing Zm), as obvious. It follows that the 
calibration of Mg2 in terms of metallicity via the comparison 
with theoretical values for SSPs is affected by a systematic 
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Mg2 

Figure 2. The difference between the average metallicity of CSPs 
and the metalHcty of SSPs having the same Mg2 index, as a func- 
tion of the index itself, for Buzzoni's 15 Gyr old models. The dif- 
ferent lines correspond to CSPs with different values for the Zm 
parameter: (0.01,0.1,0.3,0.5,1)xZq. As in Fig.l, along the lines 
Zm increases up to Z=0.1. The minimum in the line correspond- 
ing to Zm = O.O1Z0 is an artifact of the linear interpolation be- 
tween SSP models used in the computation. The effect of adopting 
different y values is shown. 



effect, which can be considerably large, if a CSP spanning a 
substantial metallicity range is present. 

The dot-dashed lines in Figure 1 show the typical range 
spanned by the observational data (e.g. Davies et al. 1987 
; WFG; CaroUo, Danziger & Buson 1993). It appears that 
galaxies with (Mg2 ,Fe52 ,Fe53 ) up to ~ (0.26,3., 2. 6) can be 
interpreted as hosting composite stellar populations with a 
metallicity distribution as in the closed box model. For these 
galaxies, higher metallic line strengths correspond to larger 
values for Zm, and because of the Mg— ct relation (Bender, 
Burstein & Faber 1993), to more massive objects. It should 
be noticed, however, that it is not possible to constrain the 
metallicity distribution from the integrated indices: the same 
value can be obtained with different metallicity ranges for 
the CSP, not for saying different distribution functions. Nev- 
ertheless, it seems that the observed metallic line strengths 
in the nuclei of the less luminous Ellipticals are consistent 
with the theoretical expectations from models for the chemi- 
cal evolution which include the occurrence of galactic winds. 

On the contrary, galaxy centers with metallic indices 
in excess than the above quoted values are not reproduced 
by these models (see also Casuso et al. 1996, Vazdekis et 
al. 1996), and appear instead to require some degree of pre- 
enrichment. The highest Mg2 are barely accounted for by 
pure SSP models, and the strongest Iron indices are not 
reproduced by the CSP model with O.IZ© < Z < 6Zq , in 
which the fraction of stars with sub— solar metallicity is less 



than 0.4. As can be seen in Figure 1, adopting a yield as 
high as 4Z0does not alter these conclusions. 

The H/9 values measured in the highest metallicity El- 
lipticals support this same picture: closed box models with 
O.OIZ0 < Z < 6Z0 have H/3 ~ 1.7, while galaxies with 
H/9 indices weaker than that are observed. Besides, the mod- 
els which account for the weakest metallic indices, with 
O.OIZ0 < Z < O.5Z0, are characterized by H/3 ~ 2, while 
galaxies with H/3 as large as 2.5 are observed. Since H/3 is 
highly sensitive to age, one can interpret the data as an ev- 
idence for a younger age of the low metallicity Ellipticals 
(Gonzalez 1993). Notice however that the high Z objects do 
require an old age. I will further discuss this point later. 

3.2 Indices versus metallicity using W models 

Figure 3 displays the models for composite stellar popula- 
tions obtained using Worthey's set of SSP models, with an 
age of 12 Gyr. The legend is the same as in Figure 1. The 
highest metallicity considered in the Z distributions is now 
~ 3Z0. The qualitative effects of the presence of a metal- 
licity spread in the stellar population are the same as al- 
ready discussed for B models, and the indications derived 
from the comparison with the observations are very similar. 
Also the quantitative effects are close to those derived us- 
ing for B models: the model with O.OIZ0 < Z < 3Zq has 
(Mg2 ,Fe52 ,Fe53 ,H/3 ) ~ (0.23,2.8,2.5,1.7) for y = 2ZQ, its 
average metallicity is [<Fe/H>] = 0.11, while a pure SSP 
with the same Mg2 has [Fe/H] =— 0.16. Compared to Buz- 
zoni's, W SSP models have a steeper dependence of Mg2 on 
[Fe/H] (see also WFG), higher Fe53 and lower H/3 over all 
the metallicity range. On the other hand for increasing 
[Fe/H], the bolometric corrections to the V and B bands 
decrease more rapidly in Worthey's models, leading to a rel- 
atively lower contribution of the high Z SSPs in equations 
( p^ to (|l3|). The two effects conspire to yield the same quan- 
titative results on the integrated indices in the CSPs. 



4 COMPARISON OF CSP MODELS WITH 
THE OBSERVATIONS 

I will now compare more in detail the predictions of the 
CSP models illustrated in the previous section to the ob- 
servations of the line strengths in the nuclei of Ellipticals. I 
concentrate on some particular aspects which can be crucial 
for understanding the modalities of galaxy formation. 

4.1 The strong Mg2 of Giant Ellipticals 

As shown in the previous subsection, a closed box model for 
the chemical enrichment in which the gas is initially at Z ~ 

fails to produce a stellar population with Mg2 as high as 
observed in the nuclei of the brightest Ellipticals. One way to 
solve the problem is assuming that the initial metallicity in 
the gas is larger than 0. To investigate this in more detail, 

1 have computed a sequence of models characterized by a 
fixed Zm = 3Z0and Z^ decreasing from Zm to O.O5Z0. The 
results are shown in Figure 4 for W 17 Gyr and B 15 Gyr 
models, where the hystogram shows the Mg2 distribution 
of the Ellipticals in the Davies et al.(1987) sample, which 
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Figure 3. The same as Figure 1, but for Worthey's 12 Gyr old SSPs. 



includes 469 objects, 
of 0.22. 



93 % of which have Mg2 in excess 



It appears that gala^cies with Mg2 higher than 0.3 (more 
than 40 % of the total sample) require Zm larger than 

O.5Z0 for both W and B models, in their nuclei. Integrated 
Mg2 indices in excess of 0.32 are obtained with Zm larger 
than ~Z0, for W, or than ~1.5Z0for B SSPs, and only 
the models computed with Worthey's SSPs account for the 
largest Mg2 . Notice that galaxies with nuclear Mg2 as large 
as 0.4 have been observed, and that the adoption of a higher 
value for the yield (dotted lines) does not solve the problem. 

The difference in the results obtained using the two sets 
of models is due to the different ages of the SSPs and to the 
different values of Mg2 at the various metallicities. Had I 



used Buzzoni's SSPs with Z in excess of 3Z0, integrated 
Mg2 larger than 0.32 would have required Z0< Z < 6Zq. 
On the other hand, using Worthey's 12 Gyr SSPs, the in- 
tegrated Mg2 for a metallicity distribution with 2Tjr..< Z < 
3Z0is ~0.33. Therefore, these two sets of SSPs require that 
the nuclei of the strongest Mg2 ellipticals host both 
old and high metallicity SSPs, with a narrow metal- 
licity dispersion. The older the population, the larger the 
allowed metallicity range, but the presence of a significant 
component of stars with Z <Z0in galaxies with Mg2 in ex- 
cess than ~ 0.32 seems unlikely, due to the strong effect it 
would have on the integrated index. 

Weiss et al. (1995) derive metallic line strengths for 
SSP models with total metallicity larger than solar, and 




Figure 4. Dependence of Mg2 on the parameter Zm: the curves 
show the loci of CSP models characterized by a maxiiimm value 
of the mctallicity Zm = 3Zq, while Zm is made to decrease. Solid 
and dotted lines correspond to y/ZQ= 2 and 3 respectively. The 
curves labelled W have been computed with Worthey's 17 Gyr 
old models; those labelled B with Buzzoni's 15 Gyr old ones. 
Also plotted is the distribution of the Mg2 observed values for 
early type galaxies from Davies et al. (1987). Objects with Mg2 in 
excess of 0.3 require Zjn greater than ~ O.SZQfor both sets of SSP 
models. 



different elemental ratios. For solar abundance ratios, their 

Mg2 indices happen to be in close agreement with those of B 
and W SSPs at Zq, but are instead systematically larger at 
supersolar metallicities. For example their 15 Gyr SSP with 
Z = 0.04 has Mg2 = 0.4, ~ 0.07 dex higher than B 15 Gyr 
model. Such high values would relieve the problem of repro- 
ducing the data for the most luminous ellipticals. However, 
it is likely that a closed box metallicity distribution with 
Zm ^ would still predict too low Mg2 values, due to the 
higher contribution of the low Z component to the optical 
light. For example, enhancing by 0.1 dex the Mg2 indices in 
W 17 Gyr old SSPs at Z >Zq, I obtain Mg2 = 0.27 for a 17 
Gyr old CSP with O.O1Z0< Z < 3Zq. The analogous exper- 
iment with B 15 Gyr models yields Mg2 = 0.30. Therefore, 
the need for a small metallicity dispersion in the nuclei of 
the strongest Mg2 Ellipticals is a robust result. 



4.2 Magnesium to Iron overabundcince 

The plot of the iron line strengths as functions of Mg2 allows 

to check the internal consistency of the indications derived 
from the different indices on the composition of the stellar 
population. In Figure 5, a selection of the models described 
in the previous section (large open symbols, see captions) is 



compared to the data from WFG and Gonzalez (1993) (small 
symbols), relative to the nuclear values of the indices. The 
solid lines show the locus described in this diagram from 
the sequence of CSP models with Zm = O.OlZeand Zm in- 
creasing up to 3Z0. It can be seen that, in these diagrams 
SSPs add as vectors: the effect of a Z distribution is that of 
shifting the model along the SSP line (WFG) , to an inter- 
mediate position between those corresponding to SSP mod- 
els with Z = Zm and Z = Zm- Besides, the line strengths 
of SSP models are stronger than those of CSP models with 
the same (average) metallicity: this can be seen in Figure 
5 comparing the position of the skeletal symbols with that 
of the corresponding polygons. Once again, this reflects the 
larger weight of the low metallicity component on the CSP 
line strengths. 

It appears that the average location of gala^des with 
Mg2 <0.3 is well reproduced by the models, better with 
Worthey's than with Buzzoni's SSPs. However, CSP mod- 
els with Zm = O.OlZebarely reach Mg2 - 0.26, Fe52 - 2.9, 
Fe53 ~ 2.7, encompassing the range occupied by the weak- 
est Mg2 objects. Galaxies with stronger metallic nuclear in- 
dices require some degree of preenrichment in their centers, 
within this class of CSP models. Notice that this is needed 
to account for both Mg2 and Iron line strenghts. 

Those objects characterized by Mg2 larger than approx- 
imately 0.3 (mostly represented in Gonzalez's sample) de- 
part from the general (Fe52 ,Fe53 ) — Mg2 relation, exhibit- 
ing lower Iron indices than the model predictions, at the 
same magnesium index. This has been interpreted as ev- 
idence for a Mg/Fe abundance ratio larger than solar in 
these (most luminous) ellipticals. According to the current 
view of how the chemical enrichment proceeds in galaxies. 
Magnesium is mostly produced by massive stars exploding 
£is Type II SNe, while a substantial fraction of the Iron is 
provided by Type la SNe. Thus, a Mg to Fe overabundance 
can be obtained either by increasing the relative number of 
Type II to Type la events (e.g. with a flatter IMF), or by 
stopping the SF process at early times, i.e. before a substan- 
tial amount of pollution by Type la events has taJsen place 
(WFG, Davies et al. 1993, Matteucci 1994). Both scenarios 
predict that all the a elements, mainly produced in massive 
stars, are overabundant with respect to Iron. It is then rea- 
sonable to assume that the [Mg/Fe] enhancement actually 
traces and a elements overabundance (or an Fe underabun- 
dance), with respect to the solar ratios. In the following 
section I will derive a simple rule which enables to estimate 
the effect of an a elements overabundance on the Mg2 and 
Iron line strenghts. 



4.3 The effect of a element enhancement on 
the metallic line strengths 

The metallic line strength of SSPs are sensitive to the metal- 
licity through two effects: one connected to the change of the 

overall shape of the isoclirone, and the other to the depen- 
dence on the metal abundance of the specific feature in the 
individual stars which populate the isochrone, described by 
the fitting functions. Both effects arc such that enhancing 
the metal abundance, the metallic features get stronger, but 
they operate in different ways. I try now to estimate how the 
shape of the isochrone on one side, and the fitting functions 
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Figure 5. Comparison between model predictions and observations of the nuclear line strengths. The data points (shown as small 
symbols) are from Worthey, Faber & Gonzalez, 1992 (crosses) and Gonzalez, 1993 (squares). The error bars quoted by the authors 
are shown in the two left panels, the smallest being relative to Gonzalez data. SSP models by Worthey and Buzzoni are displayed 
as dotted lines. The large symbols show the line strengths for a selection of CSP models: .Zni=O.OlZ0 (triangle), Zm=0.1ZQ (square), 
■Zm=0.5ZQ(penthagon) and Zni=ZQ(octagon), all having Zyi = SZQand y = 3Zq. The sequence of CSP models obtained with = 
O.OlZQand increasing Zyi is shown as a solid line. Notice that the closed box simple model fails to account for the high Mg2 and Iron line 
strenghts shown by most of the data points. The skeletal symbols show the indices of SSPs having [Fe/H] equal to the average metallicity 
of those CSP models shown as polygons with the same number of vertices. For example, the cross shows the SSP line strengths for a 
metallicity equal to the average [Fe/H] of the CSP shown as a square. 



on the other, vary in response to an a elements overabun- g hso^^^^ /c(^) Fe^2{x) dx 



dance. (^^^2 f = , (16) 



The Mg2 and Fe52 hne strengths for SSPs are given by: 



where x describes the (Log g. Log To) values along the 
isochrone, n{x) and fcix) are the number of stars and 
the contiunuum flux in the relevant wavelength bands, and 
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Mg2(x), Fe52{x) are the fitting functions, all quantities be- 
ing computed at the point x. 

For a given age, n{x),fc{x) and the fitting functions 
depend on the total metallicity Z and on the fractions = 
Xi/Z of all the different elements (Xi being the abundance 
by mass of the i-th element). The effect of changing the frac- 
tions Q (at fixed total metallicity Z) from solar ratios to a 
enhanced ratios on the shape of the isochrone is likely to be 
small. This has been shown to be valid (at low metallicities) 
for the turn-off region by Chaboyer, Sarajedini & Demarque 
(1992)and Safaris, Chieffi & Straniero (1993). Considering 
the RGB, its location essentially depends on the abundance 
of heavy elements with low ionization potentials (Renzini 
1977), which, being the main electron donors, control the 
H~ dominated opacity. The RGB temperature is then con- 
trolled by the total abundance of these elements, namely 
Mg,Si,S,Ca and Fe (Salaris et al. 1993). The a enhanced 
mixtures predicted by detailed models for the chemical evo- 
lution of galaxies (Matteucci 1992) are indeed characterized 
by a similar value of the sum of the ("i fractions of Mg,Si,S,Ca 
and Fe. For example, it varies from 0.2, at solar ratios, to 
0.18 for an a enhancement of ~ 0.4 dex. This means that a 
mixture with Z — 0.02 and solar abundance ratios has the 
same abundance (by mass) of electron donors as a mixture 
with Z — 0.022 and an average enhancement of all the a 
elements of 0.4 dex. It is then reasonable to assume that 
the shape of the isochrone is mainly controlled by the total 
metallicity, and that variations of the (i fractions have only 
a marginal effect, provided that all the a elements are en- 
hanced. Actually, compared to the solar ratios mixtures, the 
Q enhanced ones in Matteucci (1992) models have a slightly 
lower fraction of electron donors, implying slightly warmer 
RGBs and, keeping fixed all the rest, lower metallic line 
strenghts. The SSP models computed for solar ratios and 
Q enhanced mixtures by Weiss et al. (1995) indeed support 
this conclusion (compare their model 7 to 7H). 

I now turn to consider the metallicity dependence of 
the fitting functions. Buzzoni's Mg2{x) and Fe^2(x) are ex- 
pressed as the sum of two terms, one depending only on 
metallicity and the other on gravity and temperature. For 
such form, the Mg2 and Fe52 indices for SSPs can be written 
as: 

(Mg2 f = eMg([Fe/H]) + Giso (17) 

(Fe52 f = eFe([Fe/H]) + GL (18) 

in which Oivig and Gpc are exactly the dependences on [Fe/H] 
of the Mg2{x) and Feb2{x) fitting functions, while the G 
functions depend on the shape of the isochrone, and on how 
the stars distribute along it. For B models, Ofe = 1.15[Fe/H] 
and Qmcj = 0.05[Fe/H] ~ 0.05[Mg/H], since Buzzoni's fit- 
ting functions have been constructed using a sample of stars 
with likely solar abundance ratios. Asssuming that the G 
functions only depend on the total metallicity Z, and that 
the fitting functions only depend on the abundance of the 
element contributing to the feature, for B models I can write: 



(Mg2 f = 0.05 LogCMg + h(Z) 
(Fe52 f = 1.15 LogCpe + h! {Z) 



(19) 
(20) 



were all the quantities dependent on total metallicity are 
described by h, h' , and ^Mg, Cfc are the contribution to the 
total Z from Magnesium and Iron, respectively. It is now 
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Figure 6. The efTect of an a elements enhancement on the 
Mg2 and Fe52 line strenghts of W models with an age of 17 Gyr. 
The symbols are the same as in Fig. 5. The arrows show how 
the indices change for a progressively higher [Mg/Fe], according 
to the simple scaling described in the text. It appears that the 
highest Mg2 galaxies require both a larger [Mg/Fe] ratio and a 
larger total metallicity. 



easy to derive the difference between the indices for two 
SSPs with the same total metallicity (and age), but different 
a enhancements: 



A(Mg2 ) = a X LogC^Ml/Cml 
A(Fe52 ) = b X LogCg^ /^l^ 



(21) 
(22) 



with a — 0.05 and b — 1.15 for B models. 

Gorgas et al. (1993) fitting functions, used in Worthey's 
models, do not allow the separation of the metallicity depen- 
dence from that on gravity and effective temperature. Thus, 
equations (jl^) and ( |18[ ) are not strictly applicable. However, 
for bright RGB stars which dominate the Mg2 index, Gorgas 
et al. fitting functions scale according to 0Mg ~ 0.19[Fe/H] 
(WFG). Most of the contribution to the Fe52 index comes 
instead from lower luminosity RGB stars (Buzzoni 1995a), 
for which Opc ~ 1.5[Fe/H] seems a fair approximation. 
Therefore eq. (^) and ( p^ can be used to estimate the effect 
of an a overabundance on W SSP models, with a = 0.19 and 
b — 1.5, since solar elemental ratios are likely to characterize 
also the stars in the Gorgas et al. sample. 

To summarize, the line strengths of SSPs with a ele- 
ments enhancement can be scaled from those of SSPs with 
the same age and total metallicity Z (but solar elemental 
ratios) using relations (^), (^2|) provided that: 
i) the shape of the isochrone (at given age and Z), and the 
distribution of stars along it, are the same; 
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ii) the dependence of the fitting functions on the metallic- 
ity is linear in [M/H] and can be separated from the other 
dependences; 

iii) a given index depends only on the abundance of the 
element which gives rise to the considered feature. 

These three requirements are never strictly true; nev- 
ertheless within the current understanding they seem to 
be valid approximations. Chemical evolution models by 
Matteucci (1992) with an a elements enhanced mixture 
with [0/Fe] = 0.45 are characterized by ^Mg/Civig,© = 1-19 
and (^fc/Cfc,© = 0.48. It follows that, for this kind of a 
enhancement, one expects that at each metallicity point 
Mg2 increases by ~ 0.004 and Fe52 decreases by ~ 0.37 for 
B models. For W models the expected differences are larger: 
Mg2 increases by ~ 0.015 and Fe52 decreases by ~ 0.48. 
The results of the application of these scaling relations are 
shown in Figure 6, for W 17 Gyr old SSPs. One can see 
that, as long as [Mg/Fe] is fixed, the theoretical indices de- 
scribe a locus parallel to the solar ratio sequence. Therefore, 
the flattening of the Fe52 vs Mg2 relation at the high metal- 
licity end can be obtained only assuming different [Mg/Fe] 
ratios in the different galaxies. The average (flat) slope of the 
data points in the Mg2 > 0.3 domain requires an [Mg/Fe] 
increasing for increasing Mg2 , up to an overabundance of 
~ 0.4 dex. It is worth noticing that, in this simple scheme, 
Mg2 still remains an indicator of total metallicity: changing 
the [Mg/Fe] ratio at constant Z (i.e. enhancing XMg and 
decreasing Xpe accordingly), does not lead to a substantial 
increase of the Mg2 index, while the effect on the Fe52 index 
is more dramatic. As a result, the highest Mg2 objects are 
still accounted for by the highest total Z populations. Fi- 
nally, I notice that, in spite of leading to higher Mg2 indices 
for a given Z, an [Mg/Fe] = 0.4 is not sufficient to account 
for the strongest Mg2 values observed, unless the nuclei of 
these galaxies are inhabited by virtually pure SSPs of largo 
total metallicity. Supersolar CSPs barely reach Mg2 = 0.34, 
and metallicity distributions with Zm = O.5Z0don't go be- 
yond Mg2 = 0.32. Therefore, the need for a small metallicity 
dispersion in the nuclei of the most luminous ellipticals still 
remains, even assuming an overabundance which accounts 
for their relatively low Fe52 indices. 



4.4 The H/3 index 

The H/3 line strength is very sensitive to the temperature 
of turn-off stars; thus plotting H/? versus an index which is 
mainly controlled by Z may allow to estimate indipendently 
age and metallicity of a stellar system (see e.g. Gonzalez 
1993). Adopting this approach, Faber et al. (1995) suggest 
that the nuclei of elliptical galaxies form more an age se- 
quence of high Z objects, as opposed to a metallicity se- 
quence of old objects. In their comparison, the effect of a 
non zero [Mg/Fe] is taken into account plotting H/3 versus a 
newly defined index ([MgFe]), equal to the geometric mean 
of Mgb and <Fe>. This index is meant to trace better the 
total metallicity Z, but there is no guarantee that it actu- 
ally does (see also Faber et al. 1995). I prefer to use Mg2 as 
metallicity indicator, and account for the Magnesium over- 
abundance with the simple scaling given by relation (|2]]). 

Figure 7 shows the locus described by SSP and CSP 
models in the H/3 vs Mg2 plane for various ages, together 



with Gonzalez (1993) data. In order to mimic the effect of an 
a elements enhancement, a constant shift has been applied 
to the SSPs Mg2 values, which amounts to 0.015 dex for 
W and 0.004 dex for B models, coresponding to a constant 
overabundance of [Mg/Fe] = 0.4. 

The CSP models in Figure 7 appear to support Faber 
et al. (1995) conclusion, and further show that a metallicity 
spread would substantially worsen the agreement between 
the models and the observations, at all metallicities. Actu- 
ally, the dot dashed line in the left panel of Figure 7, fairly 
fitting the observations, connects the Z = 3Z0 SSPs in the 
W set of models. Since in this diagram SSPs add as vec- 
tors, an internal age spread of ~ 1-2 Gyr in the CSPs would 
hardly affect the interpretation of the data, which seem to re- 
quire younger average ages for the lower metallicity objects. 
However, since the low Mg2 galaxies are also the fainter El- 
lipticals in the local sample, a tight mass— age relation would 
be implied, with the less massive Ellipticals being (on the 
average) younger than the most massive ones. 

Basically the same conclusion holds when considering B 
models (right panel): in spite of predicting H/3 line strengths 
which are systematically higher than Worthey's, at any age 
and Mg2 , still the 15 Gyr old locus is too shallow with 
respect to the data. In essence, for old ages the model 
H/3 index is too low at low metallicities, and ages younger 
than 8 Gyr are required to fit the H/? values of the low 
Mg2 galaxies. 

The need for invoking an age difference in the galaxies of 
Gonzalez (1993) sample is related to the mild dependence of 
the H/3 line strength on the metallicity, which reflects the de- 
pendence on Z of the turn— off temperature. This is the case 
for SSP models with Red Horizontal Branches (RHB). On 
the other hand, as is well known, the HB stars in the galactic 
globular clusters become bluer with decreasing metallicity, 
although the trend is not strictly monothonic due to the sec- 
ond parameter problem (see e.g. Renzini 1977). If the aver- 
age temperature of HB stars increases for decreasing metal- 
licity, H/3 will be more sensitive to Z than estimated in the 
SSP models considered until now. As a result, the presence 
of a low metallicity tail in the stellar populations in Ellip- 
ticals could affect the CSP H/3 line strenghts appreciably, 
leading to higher equivalent widths the larger the popula- 
tion in the low Z tail. According to Buzzoni et al. (1994), 
an Intermediate Horizontal Branch (IHB) (corresponding to 
a temperature distribution peaked at Log Tc = 3.82, and 
a blue tail extending up to Log Tc = 4.05, such as e.g. in 
the globular cluster M3) leads to H/3 indices higher by ~ 
0.7 A, with respect to SSPs with RHB. Thus, in order to 
estimate the impact of this effect on the H/3 line strength I 
have computed integrated indices for CSPs using Buzzoni's 
15 Gyr old models with an artiflcially increased H/3 . The 
enhancement is taken equal to 0.5A for all metallicitites less 
than ~ O.6Z0, and linearly vanishing at Zq. The effect of 
the adoption of an IHB on the Mg2 index has instead been 
neglected, since it is estimated to lead to a decrease of only 
a few 10"^ mag (Buzzoni et al. 1994). Notice that a [Mg/Fe] 
= 0.4, which would correspond to an increase of Mg2 by ap- 
proximately the same amount for B models, has not been 
taken into account in this computation. 

The results are shown in Figure 8, where the locus de- 
scribed from CSPs with various Z^n is displayed. Obviously, 
those CSPs with Z^ in excess of ~ 0.5Zoare not affected by 




Figure 7. Compaxison between observations of H/3 and Mg2 indices in the nuclei of Ellipticals and model predictions using Worthey's 
(left panel) and Buzzoni's (right panel) SSPs. The data points (small squares) are from Gonzalez (1993) and the error bar quoted by the 

author is shown in the left panel. The dotted lines connect SSP models of constant age and different mctallicitics, each line labelled with 
its age in Gyr. The dot-dashed lines connect constant mctallicty SSPs with Z = 3ZQ(lcft panel), 1.7Zq (right panel). The solid lines are 
the loci described by single age CSPs, y = 3Zq, = O.OIZq, and Zyi increasing up to 3Z0, for the various ages. The big symbols 
indicate the location in this diagram of CSPs with fixed Zm = 3Z0and different Zm- The encoding is the same as in Figure 5, except 
for the following cases: the 3 Gyr old CSP shown as a pcnthagon has Zm = O.6Z0 (instead of 0.5Zq); the 8 Gyr old CSPs based on 
Buzzoni's models have been computed with an upper metallicity cut off of Zm ~ 1.7Z0(instead of 3Z0). This different limits are due to 
the limited Z range covered by the SSP models at these ages. A constant shift of 0.015 (for Worthey's models) and 0.004 (for Buzzoni's) 
has been applied to the theoretical Mg2 line strenghts to account for the [Mg/Fe] overabundance. 



these new prescriptions, and their corresponding line indices 
lay on the locus of RHB SSP models. On the contrary, CSPs 
with sufficiently low have substantially higher H/3 , for 
a given Mg2 , and the data in this diagram could be inter- 
preted as a seqnoiicc of old stellar populations, with increas- 
ing average metallicity. While the need for a small metallic- 
ity dispersion in the nuclei of the most luminous Ellipticals 
discussed in the previous sections leaves little room for a size- 
able contribution of a low Z component, the less luminous 
Ellipticals can host a substantial low Z, IHB component in 
their nuclei. This experiment shows that the strength of the 
H/3 line is very sensitive to the temperature distribution as- 
sumed for HB stars. Actually, just comparing Fig. 8 to the 
right panel of Fig. 7, it appears that the data can be equally 
interpreted as an age sequence, at a constant high metal- 
licity, or as a sequence at constant old age of CSPs with 
decreasing metallicity spread, corresponding to an increas- 
ing average metallicity. 



5 DISCUSSION 

The numerical experiments performed in this paper illus- 
trate how the Mg2 , Fe52 , Fe53 and H/3 line strengths are 
affected by the presence of a metallicity distribution shaped 
like the closed box model predictions. I have explored sys- 
tematically the results of changing the minimum and maxi- 
mum metallicity {Zm, Zm) characterizing the chemical pro- 
cessing: the first parameter describes the possibility of pre- 
enrichment of the gas in the closed box to different degrees; 
the second, the occurrence of galactic winds, inibiting fur- 
ther chemical processing at various levels of completion. I 
sumarize now the major results, and derive some hints on 
the stellar population inhabiting the central regions of El- 
lipticals. 

5.1 The average metallicity in the nuclei of 
Ellipticals 

Due to its major contribution to the light in the optical 
bands, the low metallicity component tends to dominate 
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Figure 8. Effect of the contribution of an IHB on the If/3 index. 
SSP models from Buzzoni with an age of 15 Gyr are shown as dot- 
ted Uncs for a RIfB and an IHB (see text). CSP model sequences 
computed with IHB SSPs are shown as different lines, along which 
Zjn increases up to 3Zq. All the CSP models have y = SZq, while 
Zin/ZQ=0. 01,0.1,0.3,0.5 for the solid, long dashed, short dashed 
and dot— dashed lines, respectively. The last CSP model of each 
sequence is shown as a big symbol. The big octagon, on the RHB 
line, shows the line strenghts for the CSP model with Zm =^0, 
= 3^0, y = 3. No shift to higher Mg2 values to account for 
the possible [Mg/Fe] enhancement has been applied here. 



the integrated indices. This imphes that the metaUic in- 
dices of the composite stellar populations are systemati- 
cally weaker than those of simple stellar populations, with 
the same metallicity. Therefore, a given value for a metallic 
lino strength, corresponds to CSPs with larger mass aver- 
aged metallicitios than SSPs. The difference between the 
two metallicitios depends on the width of the Z distribu- 
tion in the CSP model, and it can be as large as ~ 0.3 dex. 
As a consequence, a quantitative relation between metallic 
lino strength and average metallicity is subject to a substan- 
tial uncertainty: it depends on the metallicity distribution, 
which cannot be constrained, and on the SSP models, that 
may still be affected by inadequacies. Indeed, some differ- 
ences are present between the various sets of models avail- 
able in the literature (sec also Chariot, Wortlicy & Brcssan 
1996), which affect the relation between integrated indices 
and the average metallicity of a given CSP. This impacts on 
the calibration of the spectral indices of different galaxies 
in term of their metallicity, as well as on the derivation of 
abundance gradients from line strength gradients within a 
given galaxy. 



5.2 The metallicity spread in the nuclei of 
Ellipticals 

The systematic exploration of the influence of the Zm and 
Zm parameters shows that the sequence from low to high lu- 
minosity ellipticals (as far as their central stellar population 
is concerned) can be interpreted in various ways: as a se- 
quence of virtually pure SSPs of increasing metallicity; as a 
sequence of CSPs in which either the low metallicity compo- 
nent becomes less and less important, or both the minimum 
and the maximum metallicity increase. Chemical evolution 
models with galactic winds rather predict a metallicity se- 
quence among Ellipticals characterized by an increasing Zm 
for increasing mass (and luminosity) of the galaxy, with an 
important low metallicity component always present. This 
class of models can account for objects with Mg2 , Fe52 and 
Fe53 up to ~ 0.27, 3 and 2.7 respectively, while galaxies 
with higher values of these indices cannot be reproduced. 
Notice that ~ 80 % of objects in the Davies et al.(1987) 
sample has Mg2 > 0.26. Since these considerations apply 
only to the nuclear indices, classical wind models can still 
be adequate to describe the global metallicity distribution 
in Ellipticals, but a mechanism should be found to segregate 
the high metallicity component in the nuclei. 

A similar problem has been found by Brcssan, Chiosi 
& Fagotto (1994), when comparing the spectral energy dis- 
tribution of their model Ellipticals with the observations. 
These authors pointed out that the theoretical metallicity 
distribution was too heavily populated in the low metallic- 
ity range, leading to an excess light between 2000 and 4000 
A with respect to the observed spectra. They noticed that 
a much better fit was achieved if a minimum metallicity of 
Z = 0.008 was assumed, and concluded that the classical 
chemical evolution models for elliptical galaxies, like those 
for the solar neighborhood, are affected by the G—dwaf prob- 
lem, i.e. the excess of low metallicity stars predicted by the 
closed box model for the solar neighborhood. The classical 
solutions to cure the G— dwarf problem are (see Audouze & 
Tinsley 1976): 

i) infall of metal free gas, so that the SFR exhibits a max- 
imum at some intermediate epoch, when a substantial en- 
richment has been already accomplished; 

ii) prompt initial enrichment (PIE), like what explored here 
by varying the Zm parameter; 

iii) adopting a SFR enhanced in high metallicity gas, in 
which the stars are formed with a larger metallicity than 
the average Z of the interstellar medium. 

A variation of the PIE model consists in assuming that 
the first stellar generations axe formed with a conveniently 
flat IMF (Vazdekis et al. 1996) , so that they contribute met- 
als at early times, but not light at the present epoch. The fol- 
lowing generations would instead form with a normal IMF. 

All these are in principle viable solutions, and which, 
if any, of these applies to the nuclei of Ellipticals remains 
to be seen. However, I notice that the infall models predict 
for the most massive galaxies Mg2 indices not larger than 
0.28 (Bressan et al. 1996): still a low value with respect to 
the observations in the nuclei of massive Ellipticals. Solution 
iii) requires large inhomogeneities in the gas, and a variable 
IMF seems a rather ad hoc solution. 

A prompt initial enrichment for the gas in the nuclei 
of Ellipticals is easily realized by relaxing the hypothesis of 
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istantaneous complete mixing of the gas, and allowing en- 
riched gas to sink to the center. This is indeed a natural re- 
sult of dissipational galaxy formation (Larson 1975). During 
its formation process, a galaxy consists of two components: 
one dissipationloss (the newly formed stars), and one dissi- 
pative (the gas). Once formed, the stars stop participating 
to the general collapse, keeping thereafter memory of their 
energy state at formation. The gas, instead, will continue to 
flow towards the center, being progressively enriched as SN 
explosions take place. Thus, gas accumulated in the nuclear 
regions is pre-enrinched by stars which formed (and died) 
in the outer regions, and will further form stars, which will 
then be the most metal rich in the galaxy. The metal poor 
stars missing in the galaxies nuclei should be found in the 
outer regions. In other words, the different dissipation prop- 
ertities of the stars and the gas would lead to a chemical sep- 
aration within the galaxy, during its formation, no matter 
wether the protogalaxy is a monolithieally collapsing cloud, 
or if it consists of gas rich merging lumps. Interestingly, this 
out— inward formation, would also allow for peculiarities in 
the core kinematics, as observed in a substantial fraction 
of galaxies (Bender et al. 1993), since the formation of the 
nucleus would be partially decoupled from the formation of 
the rest of the galaxy. 



5.3 The elemental ratios in the nuclei of 
Ellipticals 

The Iron indices vs Mg2 plot suggests the presence of a Mag- 
nesium overabundance in the brightest Ellipticals (Mg2 > 
0.3), as pointed out by WFG comparing the data to SSP 
models. The presence of a metallieity distribution does not 
alter this conclusion, since the CSP models considered here 
describe the same locus of SSPs in the Iron vs Mg2 diagram. 
If the Magnesium overabundance is tracing a true a ele- 
ments overabundance, it is possible to estimate the enhance- 
ment quantitatively using a simple scaling of the SSP models 
constructed for solar abundance ratios. This estimate how- 
ever is very sensitive to the dependence of the Mg2 index on 
the Mg abundance, and of the Iron index on the Fc abun- 
dance. Assuming that these dependences are the same as 
the metallieity dependence of the relative fitting functions, 
for Worthey's models I have found that the brightest el- 
lipticals should be characterized by [Mg/Fe] ratios ranging 
from to 0.4 approximately, for increasing Mg2 . It may 
seem that a Magnesium overabundance would easily explain 
the strong Mg2 values without invoking high total metal- 
licities in the nuclei of the brightest Ellipticals. However, 
[Mg/Fe] >0 means Mg enhancement together with Fe deple- 
tion, with respect to the solar ratio. In the frame of an over- 
all a element overabundance, the chemical evolution models 
actually predict a lower fraction of electron donors at fixed 
total metallieity, as mentioned in Section 3.4. Correspond- 
ingly, the temperature of the RGB is increased, counterbal- 
ancing the increase of the Mg2 index conveyed by a higher 
Mg abundance. Thus, large metallicities and small metallie- 
ity dispersions are still needed to account for the data. 

According to cxirrent modelling an overabundance of 
[Mg/Fe] =0.4 implies very short formation timescales for the 
whole galaxy, since the Mg and Fe gradients, within the 
errors, are the same (Fisher, Franx & lUingworth 1995). 



How short is interesting to look at: Matteucci (1994) mod- 
els for the chemical evolution of Ellipticals predict that a 
solar [Mg/Fe] is reached already at 0.3 Gyr. An overabun- 
dance implies formation timescales shorter than this, and 
the higher [Mg/Fe], the shorter the timescale required. In- 
deed, in Matteucci (1994) inverse wind models, the galactic 
wind occurs at only ~ 0.15 Gyr in a 5 x 10^^ Mq galaxy, 
and yet the corresponding overabundance is not larger than 
0.3. The formation timescales inferred from a given [Mg/Fe] 
overabundance depend on the adopted model for the Type la 
supernovae progenitors, and accordingly can be considered 
quite uncertain. It seems however unlikely that a Mg to Fe 
overabundance could be accomplished with formation time 
scales longer than 1 Gyr. For example, Greggio (1996) 
finds that, following a burst of Star Formation, ~ 50% of 
the total Iron from the Type la SNe is released within a 
timescale ranging from 0.3 to 0.8 Gyr, for a variety of SNIa 
possible progenitors. 



5.4 The H/3 line strength as age indicator 

As for the Iron indices, the comparison of models to observa- 
tions in the H/3 vs Mg2 plot does not change if a metallieity 
spread in the stellar populations is taken into account. In 
SSPs with red horizontal branches the only way to enhance 
H/3 is by adopting a warmer turn— off, that is younger ages. 
In this case, the constraint from the metallic indices on the 
metallieity and metallieity dispersion discussed until now 
is even stronger, since younger ages make Mg2 , Fe52 and 
Fe53 weaker. Extremely large metallicities should charac- 
terize the nuclei of all Ellipticals in Gonzalez (1993) sam- 
ple, with virtually no metallieity dispersion. If however the 
temperature distribution of HB stars is wide enough, strong 
H/3 indices can be obtained in old stellar populations, due to 
the contribution of warm HB stars. In this case, the H/5 vs 
Mg2 plot would trace the different proportions of this stel- 
lar component. On the average, the effect is stronger for the 
less metallic galaxies, if they host a composite stellar popu- 
lation with larger fraction of low Z stars, and the data are 
consistent with an old age for the stars in the centers of this 
sample of Ellipticals 

The nuclei of galaxies with Mg2 in excess of ~ 0.3 are 
however likely to host only stars with metallieity larger than 
~ O.5Z0. For these galaxies, there is little room for a warm 
HB component, at least in the frame of the canonical stel- 
lar evolution, and yet their H/3 line strengths range from ~ 
1.4 to ~ 1.8 A. Adopting Worthey's 3Zq SSP models, this 
implies and age range from ~ 12 to ~ 5 Gyr, if the stel- 
lar populations in the nuclei of these galaxies are coeval. If, 
however, I consider a two component population, one very 
old and one very young, relatively high H/3 values can be 
accomplished with a small contribution to the total light 
(and even smaller to the total mass) from the young com- 
ponent. For example, a combination of a 17 Gyr plus a 1.5 
Gyr old SSP, both with 3 times solar metallieity, contribut- 
ing 80 % and 20 % of the light respectively, has H/3 ~ 1.74 
and Mg2 ~ 0.32. The contribution to the total mass of the 
young SSP would amount to only ~ 7 %. Adopting a so- 
lar metallieity for the young component, one gets H/3 ~ 1.8, 
Mg2 — 0.32 with only 10% of the light coming from the 1.5 
Gyr old population, corresponding to a ~ 3% contribution 
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to the total mass. Indeed, owing to the steep age depen- 
dence of H/3 , a small fraction of light from the young stellar 
population is sufficient to enhance the H/3 line strengh. If 
this was the case, the bulk of the population in the nuclei 
of these galaxies would be truly old, with H/3 tracing a rel- 
atively unconspicuous recent Star Formation event. 



6 CONCLUSIONS 

Using the SSP models currently available in the literature 
to construct integrated indices for composite stellar popula- 
tions with a metallicity spread I have shown that the nuclei 
of the most luminous elliptical galaxies should host stellar 
populations with: 

a) high total metallicity; 

b) a Magnesium overabundance with respect to Iron, with 

varying degrees of the [Mg/Fe] ratio. 

c) a small metallicity spread. 

Condition a) is met by processing the gas through mul- 
tiple stellar generations, and condition b) requires that this 
processing occurs within a short time scale. This inevitably 
means that during the chemical enrichment the star forma- 
tion rate was very high, implying a correspondingly large 
SNII rate. In order to proceed with the chemical processing, 
the gas had to be subject to confinement, that is it had to 
be located within a deep potential well. Condition c ) is met 
if the gas turning into stars in the nuclei of Ellipticals has 
been substantially pre— enriched, or if the maximum SFR 
was achieved at some late stage, when partial chemical pro- 
cessing had already been completed, like in infall models. 
However, for any behaviour of the SFR with time, as long 
as one considers a solfpoUuting gas mass, a low metallicity 
component in the final stellar population is unavoidable: it is 
this low metallicity component which provides the metals to 
build up the high Z stars. As a consequence, the extremely 
high Mg2 indices would rather favour the pre— enrichment 
alternative. 

These facts support the notion that the gas out of which 
the nuclei of the most luminous ellipticals formed was pro- 
duced within the galaxy itself, and was not accreted from 
the outside. Several evidences indicate that merging must 
have played a role in the formation of these galaxies, includ- 
ing the relatively shallow metal line gradients (e.g. Davies et 
al. 1993), and the peculiar kinematics in the nuclei of a sub- 
stantial fraction of galaxies (see Bender 1996). The indica- 
tions from the analysis performed in this paper suggest that 
the merging subunits should have been mostly gaseous, and 
confined within the deep potential well in which the galaxy 
itself is found. The segregation of the high Z component in 
the inner parts could result from the gas partecipating of 
the general collapse more than the stars, which could have 
formed within the merging subunits. 

As discussed in the previous section, the indications 
from the H/3 line strength are sufficiently ambiguous that 
the possibility that the bulk of the stars in the rmclci of 
the brightest Ellipticals are indeed old remains favourable. 
Therefore, the formation of these galaxies should have oc- 
curred at high rcdshifts, both the stellar component and the 
potential wells, within which the chemical processing could 
proceed up to high metallicities in short time scales. 

For the lower luminosity Ellipticals the conclusions are 



more ambiguous, as their nuclear line strengths are consis- 
tent with both wide and narrow metallicity distributions. 
However, galaxies with central Mg2 in excess of ~ 0.27 (ap- 
proximately 70% in Davies et al. sample) are not accounted 
for by closed box models with Zm ~ 0. Therefore, for the 
majority of Ellipticals some chemical separation should have 
taken place during their formation, although a substantial 
low metallicity component could be present in their nuclei. 
If this was the case, the H/3 line strength would trace the 
proportion IHB stars produced by this component of the 
composite stellar population in the different galaxies, which 
would then be old. Finally, the Mg/Fe ratio in the lower 
luminosity Ellipticals is likely to be solar, suggesting longer 
timescales for the formation of the bulk of their stars, with 
respect to the brigther Ellipticals. 

A final caveat concernes the reliability of the SSP mod- 
els used for the interpretation of the observational data. The 
differences among the various sets of models, the different 
fitting functions, and the lack in the stellar data sets of a 
fair coverage of the age, metallicity and [Mg/Fe] ratio cast 
some doubt on the use of these models to derive quantitative 
informations. The data seem to suggest that the stars in the 
rmclei of elliptical galaxies have Mg2 indices stronger than 
what can be reasonably inferred using these SSP models. Ac- 
tually Mg2 as large as 0.4 are difficult to account even with 
the oldest and most metal rich SSPs in the Worthey's and 
Buzzoni's sets of models. Although this may be a problem 
of just a few object, this fact may suggest some inadequacies 
in the SSP models. If the dependence of the Mg2 index on 
metallicity is currently underestimated at the high Z end, 
lower total metallicities, larger metallicity dispersions and 
lower [Mg/Fe] ratios would be possibly allowed for the stellar 
populations in the rmclci of the brightest ellipticals. Never- 
theless, the presence of a small metallicity dispersion in the 
nuclei of giant ellipticals and the need for old ages seem to 
be quite robust conclusions, due to the larger contribution 
to the optical light of the less metallic populations. 
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